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Abstract. Wc propose a scheme for the deterministic coherent manipulation of 
two atomic qutrits, trapped in separate cavities coupled through a short optical 
fibre or optical resonator. We study such a system in the regime of dispersive 
atom-field interactions, where the dynamics of atoms, cavities and fibre operates 
through virtual population of both the atomic excited states and photonic states 
in the cavities and fibre. We show that the resulting effective dynamics allows 
for the creation of robust qutrit entanglement, and thoroughly investigate the 
influence of imperfections and dissipation, due to atomic spontaneous emission 
and photon leakage, on the entanglement of the two qutrits state. 
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1. Introduction. 

One of the crucial ingredients in the upcoming area of quantum technologies will be 
the capability of coherently manipulating quantum systems at a distance, such that 
entanglement (i.e., quantum correlations) can be created between different nodes of a 
global quantum system. 

Entanglement is one of the most peculiar features of quantum mechanics, and 
the most distinct signature of quantum coherence. Entangled states of two or more 
particles not only play an important role for tests of quantum nonlocality [1-3], but 
also lie at the heart of quantum information processing and quantum computing [4]. 
Entangled quantum states come in many flavours, such as Bell, Einstcin-Podolsky- 
Roscn [1], Grcenberger-Horne-Zeilinger [3], or W states [5], generally depending on 
the dimensionality and tensor product structure of the Hilbert spaces involved. All 
these states have different qualities and are suitable for different roles in quantum 
information protocols [3,5]. In this context, entangled states of multiple systems 
with Hilbert spaces of dimension d (i.e. of 'qudits', with d > 2) offer their specific 
advantages over the - archetypical and most commonly considered - entangled states 
of two-dimensional systems (of 'qubits'). For instance, entangled states of two qudits 
violate local realism more strongly than entangled states of two qubits, and their 
entanglement is more resilient to noise [6]. Also, quantum cryptographic protocols 
where qubits are replaced with qudits are both more secure and faster (in that more 
information may be sent, on average, per sent particle) [7]. 

Entangled states can currently be generated in a variety of physical systems, such 
as trapped ions [8], quantum electrodynamics cavities (QED) [9], superconducting 
circuits [10], semiconductor quantum dots [11], linear optical systems [12] and impurity 
spins in solids [14]. Cavity QED [14], which concerns the interaction of atoms and 
photons inside optical cavities, provides experimentalists with a very favorable setting 
for the generation of entanglement. Atomic systems are qualified to act as qubits or, 
more generally, qudits, as appropriate internal electronic states can coherently store 
information over relatively long time scales. At the same time, in such systems photons 
are suitable for the transfer of information between distant nodes. High-finesse cavities 
can provide good insulation against the environment and can thus have long coherence 
times [15]. Two-atom Bell states, and three-particle (two atoms plus one photon) GHZ 
entangled states have been experimentally demonstrated with Rydberg atoms passing 
through a superconducting microwave cavity [16,17]. 

Schemes for the generation of entangled states of two 'qutrits' (i.e., of two three 
level quantum systems) for two atoms via a single nonresonant cavity have also been 
proposed [18]. However, in order to be used for quantum communication protocols 
[19], such an entanglement should be generated between distant atoms, like atoms 
trapped in different cavities. "Distributed" atomic entanglement requires a way to 
coherently mediate the interaction between the two atoms. One way to establish this 
interaction is through coincidence detection of photons leaking out of the cavities 
[20] : in this way, entangled states are only probabilistically generated and the success 
probability is dependent on the efficiency of photon detectors. The other possibility 
is to directly link the cavities with an optical fibre, waveguide or by a third mediating 
cavity: entangled states can be deterministically generated in such a way [21,22]. Both 
types of quantum connectivity essentially allow for the distribution of entanglement 
across a quantum network [23]. 

In this paper, we shall present a way, based on the proper choice of atomic levels' 
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structure and operating regimes, to engineer a deterministic coherent interaction 
between two qutrits embodied by atoms trapped in distant cavities, linked by a third 
optical resonator. We will then move on to study the entanglement that can be 
generated by such an interaction, as well as its resilience to imperfections and quantum 
noise. 

Before proceeding, let us first review some previous schemes for the deterministic 
generation of entangled states via such a type of connected cavities [22]. Several 
schemes have been proposed for the deterministic generation of several diverse kinds 
of entangled states [24-28], including Bell states [24], W states [28], GHZ states [26,28], 
and also qutrit entangled states [25,27]. In the schemes of Refs. [25,26], the adiabatic 
passage along dark states is employed. These schemes [25,26] are based on accurately 
tailored sequences of pulses and thus require a considerable degree of control. In other 
schemes [22,24], the Rabi oscillations of the whole system composed of the atoms, 
cavity modes, and fibre modes is utilized; the entangled states are generated through 
the exchange of excitation numbers for the atoms and photons. Hence, such schemes 
[22,24] are bound to be rather sensitive to atomic spontaneous emission and photon 
losses. In other schemes [27,28], the population of the atomic excited states can be 
effectively suppressed by virtue of dispersive atom-field interactions, but the photonic 
states in the cavities or fibre are still populated, so that the whole system is still 
sensitive to photon losses. 

The scheme we propose here is different from all such previous schemes [22,24- 
28], even from those adopting adiabatic passage [25,26]. The scheme is inspired by 
a previous idea for virtual-photon-induced phase gates between two distant atoms 
[29]. In the scheme, the entanglement is generated through virtual population of not 
only the atomic excited states but also the photonic states in the cavities and fibre. 
Therefore, our scheme will turn out to be well shielded from both atomic spontaneous 
emission and photon losses. 

The paper is organized as follows. In Sec. [2j we specify our conditions on the 
physical parameters and derive the effective Hamiltonian for the system. In Sec. |3l we 
discuss the generation of two qutrit entangled state via the effective Hamiltonian 
and study the reliability of the entangled state in the presence of mismatches in 
the system's parameters. In Sec. [4j we discuss the influence of atomic spontaneous 
emission and photon leakage on the entangled state. Sec. [5] contains some concluding 
remarks. 

2. The model 

The schematic of our setup is shown in Fig. 1. Two distant atoms are individually 
trapped in two double- mode cavities (A and B), which arc connected by a third optical 
resonator, as shown in Fig. 1 (a). The linking resonator can be either a third cavity 
coupling the two distant cavities (like in a photonic crystal), or a 'short' (in a sense 
which will be specified shortly) optical fibre. For simplicity, we will henceforth refer 
to the linking resonator as to the "fibre" . 

The coupling of the fibre modes to the modes of the cavities in Schrodinger interac- 
tion picture may be modeled by the Hamiltonian Hj = Yl^Li ^2k=L r ^-n,kb\ L fhn.k + 
v n ,k{b n ,kW A k + {—l) n e itpt o) B k )] + H.c.} (H = 1 is used throughout this paper), where 
A nj fc is the frequency difference of the nth polarised fibre mode and the cavity mode 
with the corresponding polarisation [L and R denote, respectively, er^-circular and 
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Figure 1. Setup and atoms' levels' configuration for realising qutrit 
entanglement, (a) Two atoms are trapped in the double-mode cavities A and 
B, respectively; the cavities are coupled by an optical fibre. (b) Possible 
implementation with 87 Rb atoms, showing the involved atomic transitions for 
each atom. 

er~-circular polarisation), b n ,k and a^.fc ( a s,fc) & re the annihilation operators for the 
polarised modes of the fibre and of cavity A (B), v n ^ is the corresponding coupling 
strength, and the phase cpf is due to the propagation of the field through the fibre 
of length L: ip = 2itujL/c [22]. In the short fibre limit 2Lv/(2ttc) < 1 [22], where 
77 is the decay rate of the cavity fields into a continuum of the fibre modes, only the 
resonant modes &l and bu of the fibre are excited and coupled to the cavity modes. 
In this case, the interaction Hamiltonian H C j describing the cavity-fibre coupling can 
be rewritten as [22,27] 

H? = Mh(a\ :k + e^f^ Bk ) + H.c.]. (1) 

k=L,R 

In this paper, the state of the photon modes for cavity A (B) or the fibre is taken to 
be \ii') s (s = Ci, C2, and fib), with i (i r ) denoting i eft- {if cr~-) photons. 

The atoms have three excited states (le^), |e ), and \e R )) and three ground states 
(I<?l)i I <?o) : an d 1.9/?)), which could be the Zeeman sublevels of alkali atoms in the 
excited- and ground-state manifold, respectively. To fix ideas and portrait a case 
of practical interest, we consider here a possible implementation with s7 Rb, whose 
relevant atomic levels are shown in Fig. 1. (b). We only illustrate the involved state 
transitions by starting from the initial state \go)A \9o)b f° r the atoms. Each atom is 
assumed to be coupled (off-resonantly) to an external 7r-polariscd classical field and 
both a'- and cr~-polarised photon modes of the local cavity. 

We first describe the involved transitions for each atom in its local cavity. In cavity 
A, the transitions \go) — * |eo) and |eo) — > \gi) (!(?#)) are coupled to the 7r-polarised 
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classical field and the er^-circular (er~ -circular) polarised cavity mode, respectively. 
In cavity B, the transitions \g ) — > \e L ) {\e R )) and \e k ) — > \g k ) arc coupled to the 
er^-circular ((T~-circular) polarised cavity mode and the 7r-polarised classical field, 
respectively. It should be noted that, for the selected frequencies of the classical and 
cavity fields, additional transitions cannot be induced due to the large difference of 
the energy levels between the F = 1 and F — 2 states of the ground manifold h 2 S\/2- 
In interaction picture, the Hamiltonian describing the interaction of the atoms with 
the cavity and classical fields can then be written as 

Hf d = E IdA^ke 1 ^^ \e ) A (g k \ + fUe^-W^ ^ ^ 

k=L,R 

+ 9B,ka B ,ke iAB ^ mt \e k ) B (g \ + £l B e iAa "»"> t e i + B \e k ) B (g k \ + H.c], (2) 

where A I;JZ (x = A, B; y = eQ,eL,e-R\ z = ga,gL,<jR) denotes the energy difference 
between the fields and the corresponding atomic transition \y) «-> \z) in cavity x; g x<k 
is the coupling strength of the atom with the polarised photon mode in cavity x and 
satisfies g Xfk = <?oCm,m' (with go and C mjTO ' being the atom-cavity coupling constant 
and Clebsch-Cordan coefficient, respectively); fi x and <fi x are one-half Rabi frequency 
and phase of the classical field, and H.c denotes Hermitian conjugate. 

Under the condition of larger detuning, i.e., A x ,yz 3> 9x,k, the probability 
that the excited atomic states are populated is virtual, then the Hamiltonian (2) is 
reduced to an effective one that involves only the Stark shifts induced respectively by 
the classical and cavity fields for the three ground states, and Raman transitions 
l<?o) — ► \g k ) {k = L, R) induced collectively by the classical and cavity fields 
(See Appendix A). Furthermore, to avoid excitation of real photonic states in the 

cavities and fiber, we first set ii\ = — = A ° B — , u 2 = A g ' 4,fc — = , 9s ' k — , 

\ — gA.k^A I 1 | 1 \ 9B,k^B I 1 | 1 ^ A — A A 

A " 2 ^AlJ^ + ^I^) - -^iA^^ + A^:)' A = &A,e ogk -&A,e ogo ~ 

&-B,e k g — ^B.e k g k , 4> 'A = 4>b an d v = v k , and satisfy the condition \[2v 1 |A — V2v\, 
A + y/2v, and A > f, §. In this case, the energy exchange between the atoms 
and the photonic modes of the cavities and fiber is also virtual (see Appendix A). 
Suppose all the modes of the cavities and fiber are initially in the vacuum state, i.e., 
|00) Ci \00)j: ib 1 00) C2 . Thus all these modes will remain in the vacuum state during the 
evolution. Therefore, the global effective Hamiltonian reads [29] 

K= E v(\9o) A (go\ + \g k ) B (g k \)-x(e- llpf S Ak S B:k + H.c), (3) 

k =L,R 

where 

A2 r/ 1 1 2 n 

" = ^ + ^x^ + Awr, + ^ (4) 



A 2 , 1 1 2 



4 A - V2v A + V2is A ' 



(5) 



S A k — \go) A (g k \, and S Bk = \go) B (g k \- The Hamiltonian (3) allows for the global 
coherent manipulation of the atomic states. We will show this in detail by studying 
the generation of qutrit entanglement between the two distant atoms. 
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3. Generation of qutrit entanglement 

We now show that the effective Hamiltonian (3) allows one to generate a qutrit-qutrit 
entangled state between two atoms A and B. Initially, the two cavities and the fibre are 
in the vacuum state while the two atoms are initialised in \ip AB (0)} = \go) A \go) B (this 
can be achieved by optical pumping with two classical laser fields, one at resonance 
with the transition from F = 2 to F' = 2, the other one coupling to the transition 
from F = 1 to F' = 2 [30]). From the effective Hamiltonian (3), we immediately 
obtain the temporal evolution of the two atoms as follows: 

\^AB(t)) = e -^[cos(V2 X t) \g ) A \g ) B 

+ -j=e~^ S m(V2 X t)(\g L ) A \g L ) B + \g R ) A \g R ) B )]. (6) 

Setting yt = a ' ctan (jj0+ m7r ( m = 0,1,2,...), we get a qutrit-qutrit maximally 
entangled state (in the sense that the local Von Neumann entropy is maximal) 

= & -^-\\9o) A \go) B + ie- i{p t{\g L ) A \g L ) B + \g R ) A \g R ) B )]. (7) 

The scheme is deterministic and viable for a rather wide range of system parameters. 
It should be noted that the process for the generation of the qutrit-qutrit entangled 
state in the present scheme is at variance with previous strategies adopting dispersive 
interactions [27], because here the occupation of both the atomic excited states and 
of the photonic states in the cavities and fibre are negligible. Moreover, the external 
control in the present scheme is less demanding, as the preparatory step to put the 
atom in a specific superposition of two ground states [25,27] or the local manipulation 
of one of the atoms (by yet another classical field) during the temporal evolution of 
the whole system [25], are not required here. 

In the above analysis, exact knowledge of the system parameters is assumed. 
But, in general, there could be various errors in the parameters due to the imperfect 
characterisation of the system. The potential errors include: 

• the mismatch of the coupling rate g x .k and Q x (x — A, B\ k = L.R.) for the 
atoms with the local cavity and classical fields, as g x ^ and £l x are dependent on 
the atomic position and might fluctuate; 

• the mismatch of the phase <p x due to noise in the phases of the classical fields; 

• the mismatch of the detuning A x ^ yz (x = A, B; y = eo, e^, e R ; and z = go, gL, gn) 
between atoms and fields due to possibly imprecise control; 

• the mismatch of the coupling rate for the cavity and fibre modes, as Vk is 
decided by the manufacture technology and might be imprecise; 

• timing errors, due to the finite switching rates of the interactions and the limited 
precisions of the interaction times (t A and t B will be set for each atom interacting 
with the local fields.); 

• polarisation errors due to unstable magnetic fields, which lead to mismatches 
in two parameters related to polarisation, such as g x ,L(R)> v l(r) an d A X ,L(R) 
(A x L ( B <j here denotes the detuning for it — or tt — o~ Raman channel). 

In order to check out how the mentioned errors influence the generation of the 
entangled state, we define the following fidelity as a measure of the reliability of the 
qutrit-qutrit maximally entangled state: 

F = ^\ D B\Tr Cl j :C Mt)}\^ A D B ) , (8) 
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Figure 2. The fidelity of the qutrit-qutrit entangled state versus kinds of errors 

(all the parameters plotted are dimensionless). (a) F vs 3A < L W and Sb,l{r) 

9A,L(R) 9b, hr) 

( b ) F VS S3A (B),L *9A(B),R (c) p vs SOa and SOb {d) F vs an( J % 

9A(B),i 9A(B),R "A v ' 0A *B 

(e) F vs and (f) Fvs^ and (g) F vs and . and 

( h ) F vs »fA(B),L and '^(B),^ 



where p(t) is the state of the entire system at arbitrary time (governed by Eq. (9), 
where neither decoherence nor errors are accounted for), and Tr cl j iC2 denotes the 
partial trace over the field degrees of freedom. 

We first assume "perfect" interactions, considering the case g x> k = g, Q x = ^ = g, 

^A,e g = A B ,e k g k = 20g, A A ,e g k = ^B,e k g = 21g, <j> A = 4>B , and V k = V = y^g) 

as a reference. Under such conditions maximal qutrit-qutrit entanglement is obtained 
at the reference time t = To (in the notation of the previous section, only the case 
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m = is considered, i.e., xTq = 0.6755.). We then set the errors involved in the 
parameters g x>k , Cl x , A x ^ yz , v k , (f) k and t x to be 6g x>k , 6ft x , SA X)VZ , 5v k , 6<f> k and 
St x , respectively. In Fig. 2, the fidelity is plotted versus the different kinds of errors. 
Notice that these fidelity plots display a number of symmetries. The mirror symmetries 
about the line bisecting the axes trivially reflect the choices of the error parameters 
and the symmetry of the system under exchange of the two atoms and cavities. Some 
of the symmetries arc instead more interesting: for instance, errors in the detunings 
induce additional phases and have clearly oscillatory effects (g), errors in the cavity- 
fibre coupling strengths induce different Stark shifts (and have hence different effects) 
depending on their signs (f), while errors in the atom-light coupling strengths in 
the two different polarisations have, more intriguingly, approximately 'rotationally 
symmetric' effects (b). Let us now quantitatively discuss the influence of the various 
kinds of errors on the fidelity with the entangled 'reference state'. 

It can be seen from Fig. 2 (a), (6), (e) and (/) that the fidelity F is very robust 
against errors in the parameters g x , k , v k and t x . A deviation |<5g Xi fc| ~ 10%g x ^ k , 
\8v x \ ~ 10%f x , or \St x \ ~ 10%To will cause only a reduction smaller than 10~ 2 in the 
fidelity. 

From Fig. 2 (c), it is apparent that the fidelity is, on the other hand, very sensitive 
to imperfections in n x , mainly dependent on the Stark shifts induced by the classical 
fields [see Eq. (3)]. However, the influences of imperfect fl x through such Stark shifts 
can be eliminated as one can apply a second classical field to produce offsetting ac- 
Stark shifts on both atoms [27]. If this is done, then the effect of the errors \S£l x \ on 
Q x will be analogous to the effect of the deviation |&7a;,k| in g Xyk , which have already 
been shown to be very slight. Thus, this simple countermeasure would make the 
entanglement fidelity robust also against possible errors in Cl x . 

When deriving the effective Hamiltonian (3), we set the condition </>a = <t>B to 
eliminate the phases of the classical fields in Eq. (3) (see Appendix A). However, 
such phases are affected by noise and fluctuations, and might be slightly different 
in practical instances. Nevertheless, the fidelity is only marginally degraded by the 
possible errors in the parameter 4> x : a deviation \8(j) k \ = 3%(p k will cause only 2% in 
the reduction of the fidelity. In practice, Using only one classical field to illuminate 
both atoms that are distributed in the two cavities would help in keeping the phase 
fluctuations under control [31]. 

Let us remind the reader that in our scheme the cavity and the classical fields are 
detuned from the corresponding atomic transitions by specific values. It can be seen 
from Fig. 2 (g) that the fidelity is highly dependent on the parameter A x ^ yz : a small 
deviation in A x , yz leads to large oscillation in the fidelity. This is mainly due to the 
detuning-dependent Stark shifts induced by the cavity and classical fields, besides the 
possible occurrence of a phase (A = A A . eagk - A A .e aga - A B>ekgo + A B , ehgk ) in the 
exponential factor e tAt in Eq. (3). Though this requirement is strict, it is not a major 
problem with the currently developed laser technology in cavity QED experiments: the 
necessary stabilisation of the fields' frequencies can be achieved by means of acousto- 
optic modulators [32]. In Fig. 2 (h), we check the stability of the fidelity versus the 
detunings for the two polarisation (tt — cr or 7r — a~) Raman channels. Ideally, the 
conditions Aa.l=Aa,r and Ab i l=Ab i r are required. But, in real experiments, this 
requirement may not be perfectly satisfied, due to the fact that the magnetic field can 
break the degeneracy between the atomic ground states. Our investigation, reported 
in Fig. 2 (h) shows that the fidelity is only slightly affected by errors in the detunings 
for both polarisation channels: the fidelity with the entangled state of reference will 
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Figure 3. The fidelity of the qutrit-qutrit entangled state versus the 
dimensionless parameters \t, K/g, "y/g or &/ g. (a) F vs \t and K/g (re = (3 = 7); 
(b) F vs K/g and -y/g ((3 = 10~ 3 g); (c) F vs -y/g and j3/g (re = lCT 3 g); and (d) 
F vs K/g and /3/g (y/g = lO" 3 ^. 

still be larger than 0.96 even when a deviation tfA^.s)^ ~ 3%AA(B).k occurs. 

4. Influence of spontaneous emission and photon leakage 

In all the above arguments, we have assumed the entire system is ideally isolated from 
the environment, and have not considered any dissipation. In this section, we take 
into account the dissipation due to atomic spontaneous emission and photon leakage 
from the cavities and fibre. The master equation for the density matrix of the entire 
system can be expressed as 

f>= -i[H fulh p] + - [ {2a X}k pal tk -al k a X)k p-pal k a Xik ) 

k=L,R x=A,B 

+ ^(2b k pb{-blb k p-pblb k )] 

+ 5 Yl Yl (2A x , aP Al a - Al a A x , aP - P Al a A x , a ), (9) 

x=A,B tr=L,R,n 

where A Xt „ = Y, y , z \v) x (2/5 1(7 I z ) x ( z \ (v = 9l, go, 9r; z = e L , e , e R .) is the atomic 
lowering operator, with x (y; la \ z) x being the Clebsch-Gordan coefficient (i.e., C mj7Ti ') 
for the dipole transition |e) <-> \g) with polarisation a = L, R, tt; 7, j3 and k stand, 
respectively, for the spontaneous emission rate and for the fibre and cavity decay rates 
(assumed for simplicity to be equal for the two cavities and for the two polarised 
modes). The contribution of the thermal photons have been neglected, as is possible 
at optical frequencies. 

The master equation (9) has been numerically solved in the subspace I G {r f u u , 
\9l) a \9o) b |00) Ci |00) /i6 |00) C2 , \g R ) A \g ) B |00) Ci |00} /i6 |00) C J. In Fig. 3 (a), the 
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fidelity of the maximal qutrit-qutrit entangled state is plotted versus the dimensionlcss 
parameters \t an d ft/g (k — (3 — 7 is set). In Fig. 3 (b), (c) and (d) the 
fidelity is plotted versus each pair of the three dimensionlcss parameters n/g, (3/g 
and 7/5 (the remaining one is set to be 10~ 3 g). In the calculations, we still set 
fl A = £l B = g A<k = g B<k = g, A AiBogo = A B , ekflfc = 20g, A A>eogh = A B ,e kgo = 21g, 
and Vk = v2g. 

From Fig. 3 (a), we note that the fidelity is almost unaffected by the three decay 
rates k, (3 and 7 when k = (3 = 7 = 10 g. Even when k = (3 = 7 = 10~ 2 g, the fidelity 
is close to 0.97, which is much larger than the one (< 0.87) obtained in Ref. [27]. This 
improvement is of course due to the suppression of the excited states' population of 
the fields, as well as of the atoms. From Fig. 3 (6), (c) and (d), it can be seen that 
a decay rate of 10~ 2 <?, for either n, (3 or 7 alone (with the other two parameters set 
to zero) leads to a fidelity larger than 0.98. Note that the previous scheme through 
the adiabatic passage [25], a decay rate n = 10~ 2 <7 alone degraded the fidelity down 
to F = 0.95. 

5. Conclusion 

In summary, we have proposed a scheme of atomic levels (with an explicit possible 
realisation in Zeeman sublevels of alkali atoms), where qutrit quantum information can 
be stored in three ground states and, most importantly, manipulated globally between 
distant nodes through the virtual excitation of excited atomic levels and mediating 
bosonic fields like, typically, light. 

Our scheme is different from any previously proposed ones in that this choice of 
atomic levels allows for the whole coherent evolution of the global system - involving 
the two atoms and the linking bosonic modes - to be driven by the virtual excitation 
of both the atomic excited levels and the intervening fields. 

This feature renders our scheme remarkably more robust than any other 
previously proposed in the face of decoherence, whose main sources in these settings 
are photon loss and spontaneous emission from excited levels. Also, our scheme - 
being based, essentially, on the proper choice of atomic levels and operating regimes 
- requires very modest control and proves to be rather resilient against experimental 
imperfections as well. All these qualitative remarks have been substantiated in this 
work by a very thorough quantitative analysis of such unwanted effects. 

Clearly, a price has to be paid for improved robustness: the use of exclusively 
virtual excitations makes these coherent manipulations very slow if compared to 
schemes adopting resonant couplings [33]. Ultimately, the choice between faster, 
resonant schemes and more robust, virtual ones should depend on the use one intends 
to make of them. Of course, speed would be paramount in applications directly 
related to quantum computation. However, as pointed out in the introduction, qutrit 
systems are mainly interesting for quantum communication purposes, where the most 
delicate task to accomplish is precisely the robust distribution of entanglement between 
distant nodes of a network, and speed is not as crucial: the present, fully virtual 
scheme would respond precisely to this need. In this perspective, our study shows that 
atomic systems hold considerable promise for the encoding and coherent, distributed 
manipulation of multidimensional quantum alphabet for quantum communication 
purposes. 
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Appendix A: Derivation of the effective Hamiltonian 

Under the condition of large detuning, i.e. for A XtVZ ^> g x ,k, Qx, and as long as 
the atoms are initialised in the ground states, the probability that the atomic excited 
states are populated is virtual. Thus, the atomic excited states are negligible during 
the time evolution of the entire system. In this case, we can adopt the time-averaging 
method [34] to obtain an effective Hamiltonian as follows [35]:. 



Hf*= -iHf d (t) J Hf cl (t')dt' 



im A \9q)a (sol + m B \9k) B (gk\ 



k=L.R 

+ Va A , k a A ,k a A,k \9k) A (9k\ + ^a B , k ^B,k a B,k \9o) B (9o\ 

+ K A:k n A (aA, k e lAA < kt e-^ A \g ) A (g k \ + H.c) 

+ Xa^ k n B (a B Me lAB - kt e~ 1 ^ \g k ) B (g \+H.c)}, (10) 

where M tu = sij^, Mn B = s^J^ ^ A , k = Va B , k = sgj^, K A , k n A = 

9A,k&A I 1 I 1 \ \ gB.fc^B I 1 , 1 \ A A 

—2— \ + A A , eo9 J' <Wfi B ,* - 2 lSi^ + A B . efc9 J> A ^ fe - A ^o 9fc - 

Aa, EoSo , A B , fe = A B ,e h g a - A Biefc9fc . For H£ ci in Eq. (10), the first (second) and third 
(fourth) terms describe the Stark shifts for the states |<7o) (|<7fc)) and \g k ) (\go)) 01 the 
atom in cavity A (B), induced by the classical and cavity fields, respectively; the fifth 
(sixth) term describes the Raman coupling between the states \go) and \g k ) for the 
atom in cavity A (B), induced collectively by the classical and cavity fields. 

Hence, the effective Hamiltonian of the entire system is given by H e = Hj + H^ cl . 
Let us now introduce three normal modes c k and c± k by applying the canonical 
transformation c k = -^(a^.fc - e~ lipi a B , k ) and c ±k = \{a A . k + e~ i<p faB,k ± V2b k ) 

[22]. Then, we switch to a rotating frame by the unitary transformation R = e~ lS ' ' 
[29], i.e., H' e = R)H e R - iR^, and obtain 



K = {^Qa \9o) a (5o I + m B \9k) B (9k\ 



k=L,R 



+ ^^(c+feC+fe + c^ fc c-fc + 2c\c k ) \g k ) A (g k \ 
+ -^(c+fcC+fe + ci fc c_ fe + 2c\c k ) \g } B {go\ 

+ ^( C \ k c^^ + V2c\ k c k e^t + Sl c \_ kCke -iV2»*t + Hx) (fffc | 
+ ^L[ c \ kC _ k e^^ V2c\ k c k e^ - V2cl k c k e~^ + H.c) \g ) B (go\ 
^-[(c+ke-^* + c^e 1 ^ + V2c k y A -^e-^ \g ) A (g k \ + H.c] 



[{c+ke-^ 1 + c. k e^ 1 - V2c k )e lA -^e-^-^^ \g k ) B (g Q \ + H.c]}. 



^a B t.n B 



2 

For simplicity, we now set ^ = (j,q a = /j,n B , /i 2 = /i aA>fc = Ma Bjfc , A = X aA , k U A = 
\a Bk n B , A = A A ,k = A Bk , 4>a = 4>b and v = v k . Under the condition \f2v 
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|A - y/2v\, A + y/2v, and A > f, §, the energy exchange between the bosonic 
modes and the atoms as well as between the bosonic modes themselves is virtual. The 
virtual excitation of the bosonic modes leads to the Stark shifts and coupling between 
the atoms. Then H' e reduces to [29] 



K = M\9o)a (5ol + \9k) B (9k\) 



k=L,R 



+ ^( c +fe c + fe + c -k c -k + 2c l c k)(\9k) A (9k\ + \go) B {go\) 
-( c lk c +k - cl k c-k)(\gk) A (gk\ + \go) B (go\) 2 



32V2V 

-^-{c\ k c +k - ci k C- k )(\g k ) A (g k \ - \g Q ) B (g \) 



1 12 



c\ k c +k + cl-fcC-fc + -rclc k ){\g k ) A (g k \ + \g ) B (g \)} 



A - V2v A + V2v ~ k A 

- 4"( X —F= + i)(e~ lvl S\ k S R k + H.c)}, (11) 

4 V A-y/2v A + V2v A A Ax B ' k ;/ ' v ; 

with S A k = I <7o) a (9k\ and S B k = \go) B (g k \- The quantum-number operators c^ +k c+ k , 
c lfe c -fci CfcCfc for the bosonic modes are conserved quantities during the interaction as 
all of them commute with the Hamiltonian H". Suppose all the modes of the cavities 
and fibre are initially in the vacuum state, i.e., |00) Ci |00)y ib 1 00) Ca . Hence, all the 
bosonic modes c+fc, C- k and c k will remain in the vacuum state during the evolution. 
Finally, the global effective Hamiltonian H'J reads 

K= E v(\go) A (g \ + \gk) B (gk\)~x(e- lipf sl k S B , k + H.c), (12) 

k=L.R 

where 

V = l l i+V / , (13) 

"'4l<A^ + ZT^ + !>- <14) 

and 

^^-x^-a-ttt^^- (15) 

It should be noted that we have neglected the term rf \go) A (go | to maintain the 
symmetry in the effective Hamiltonian (12). In practice, this term can be compensated 
by an additional ac-Stark shift for the state \go) of atom A [27]. 



Appendix B: Validity of the effective dynamics 



We now turn back to the full Hamiltonian of the system and check how accurate is the 
description of the system through the effective Hamiltonian (3). We take the energy 
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level \F = 2) of 5 2 Si/ 2 to be the zero energy reference point, and write down the full 
Hamiltonian for the entire system as follows: 

Hf u u = Y (ajf,kblb k + Lu aAk a Ak a A , k + u aB ,k a B,k a B,k 

k=L,R 

UA, 9k \gk) A (g k \+ujB,e k \ e k ) B (e k | ) + u A , ea leo)^ (e | + ub, 9o |5o) B (go\ 
+ Y l9A,ka A ,k \e ) A (g k \ + n Ae -^^-M \ eQ ) A {go \ 

k=L,R 

■ II /•;./, " /;./, \e k ) B (ffo| + n B e-< Ua B*-M \ e ^ g (g k \ + H.c\ 
+ Y Mbkia^ + e^fa^ + H.c], (16) 

k=L,R 

where cj/^ and u aA k {u> aB k ) are the energy levels for the polarised photons in the fibre 
and cavity A (B), respectively, wq x (a; = A, B) denotes the energy for the 7r-polarised 
classical field Q x , and u> x m (m = g , gk, eo, ej.) is the energy level for the atomic state 
\m) x . 

Taking the initial state | -0(0)) = \go) A |<?o) b 1 00) Ci |00)y ib |00) C2 and considering 
all possible states of the system in evolution, we express the state of the system at 
time t as \ipf u ii(t)) — Ci{t) l^i) ( c i(t) being time-dependent amplitudes) within the 
subspace T f u u spanned by the vectors {\<pi) , • ■ ■ , \4>i) , ■ ■ ■ , |</'i2)}: 

Ffuii = {(\go) A \go) B , \e ) A \go} B , \9l) a \gi) B , \gn) A |S-R> B > \9l) a \e L } B , \9r) a \zr)b) 
®|00) CI |00) /4b |00) C2 ,| 5i ) A |. 9o ) B 

® (|10) C1 |00) /ffi |00) C2 , |00) C1 |10) /46 |00) C2 , |00) ci |00> /46 |10) C2 ), \g R ) A \g ) B 
® (|01) C1 |00) /16 |oo) C2 , |00) ci \oi) fib |00) C2 , |00) ci \Q0) fib |01) C2 )} . 
The occupation probability for each state vector |^) during the evolution is P,(t) = 
\ci(t)\ , and satisfies J2i=i ^»(*) = 1- Thus the occupation probability of the atomic 
excited states and the photonic states arc P e (t) = J2i=2 5 6-^ anc ^ Pp(t) = Si=7-^> 
respectively. The validity of the effective Hamiltonian implies that both the occupation 
probability P e {t) and P p (t) should be small enough thus they can be negligible during 
the time evolution of the entire system. We focus here on quantum state transfer, 
i.e., on the variation of the occupation probability P\(t) and Pt ra {t) = Sj=3 4-Pij 
for the numerical verification of the effective dynamics, which is portrayed in Fig. 4. 
Fig. 4 (al) and (a2) are obtained through the solution of the Schrodinger equation 

i^u!M. = H'J \ipeff(t)) in the subspace T e € {{fa), \<t> 3 ), \4>i)}- The two figures 
display perfect Rabi oscillations, which indicates ideal state transfer between the states 

|^) and 1/V2(|0s> + |&)). 

Fig. 4 (61), (cl), (dl), (62), (c2) and (d2) are obtained through the solution of 

the Schrodinger equation i ^rtl^iQl = Hf u u \ipf u ll(t)) in the subspace Tf u u. It can be 
seen from Fig. 4 (61) and (62) that state transfer via the effective Hamiltonian (3) is 
almost perfect (for Fig. 4 (62), this is especially apparent.), indicating that numerical 
results obtained from the effective and full Hamiltonians would be equivalent. Fig. 4 
(cl) and (c2) plot the variation of the occupation probability of the atomic excited 
states, while Fig. 4 (dl) and (d2) plot the variation of the occupation probability of 
the photonic states in the cavities and fibre. It is apparent that both P e and P p are 
very small during the evolution of the entire system. The analysis above verifies the 
validity of the effective Hamiltonian (3). 
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Figure 4. The occupation probability Pi(t), Ptra(t), Pe(t) and P p (t) versus 
the dimcnsionlcss parameter gt, respectively, Qa = = gA,k = 9B, k = 9i 
AA,e M o = &B, £k g k = Ai, AA, eo9fc = A B ,e k g = A 2 and 1/ = v^g- (al) ~ (dl) 
Ai = 10g, A 2 = llg; (a2) ~ (d2) Ai = 20g, A 2 = 21g. 



Let us now review the physical conditions given in the description of our system, 
to reveal some insight about and relationships between certain dynamical parameters 
when the effective Hamiltonian (3) is valid. In order for H" to hold, we required 
A, > g, n as well as y/2v, |A± y/2v\, A > ^ (A < A h i = 1,2). A, is 
the dominant factor, because the occupation probability of the atomic excited states 
(P e ) and the photonic states (P p ) are inversely proportional to Af, given all other 
parameters are pre-set. This is proved to be true in Fig. 4 (cl), (dl), (c2), and (d2), 
as the average occupation probability P e and P p shown in Fig. 4 (cl) and (dl) are 
about four times that in Fig. 4 (c2) and (c?2). In other words, the difference between 
the effective Hamiltonian (3) and the full Hamiltonian (16) decreases with increasing 
A^ This can also explain the phenomena for the different deviation from the perfect 
state transfer, which are more noticeable in Fig. 4 (61) than in Fig. 4 (62). 

Let us stress once again that the obtained effective Hamiltonian (3) is indeed 
valid as the occupations of the atomic excited states and the photonic states have 
been showed to be strongly suppressed. 
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